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ABSTRACT 

We analyze scale dependence of redshift space bias b and (3 = Vt^-;^ /h in the context 
of the halo model. We show that linear bias is a good approximation only on large 
scales, for k < O.l/iMpc"^. On intermediate scales the virial motions of galaxies cause a 
suppression of the power spectrum relative to the linear one and the suppression differs 
from the same effect in dark matter. This can potentially mimic the effect of massive 
neutrinos and the degeneracy can only be broken if power spectrum is measured for 
k ^ 0.1hMpc~^. Different methods to determine /3 converge for k < 0.1hMpc~^, 
but give drastically different results on smaller scales, which explains some of the 
trends observed in the real data. We also asses the level of stochasticity by calculating 
the cross-correlation coefficient between the reconstructed velocity field divergence 
and the galaxies and show that the two fields decorrelate for k > 0.1hMpc~^. Most 
problematic are galaxies predominantly found in groups and clusters, such as bright, 
red or elliptical galaxies, where we find poor convergence to a constant bias or /? even 
on large scales. 



1 INTRODUCTION 



Determination of the power spectrum of mass fluctuations 
is one of the main goals of existing and u pcoming galaxy 



surve ys. Current state of the art is PSCz (Saunders et al 



200C), which has a near spherical geometr^r and consists of 



about ll.'^nnn mpasiTrprl galayy rpHshiftt; TTprnmrng fiurvpyK 



such as 2 degree Field survey (2dF)[j and Sloan Digital Sky 
Survey (SDSS) Fl will measure redshifts of up to a million 
galaxies. 3-dimensional mass power spectrum is sensitive to 
a number of cosmological parameters, such as the matter 
and baryon density, shape and amplitude of initial fluctu- 
ations and the Hubble constant. This sensitivity is further 
improved if additional information from cos mic microwave 



background (CM B) anisotropics is included (Eisenstein, Hu 
& Tegmark 1999). Mass power spectrum is particulary im- 
portant for determination of neutrino mass. Massive neutri- 
nos have only a minor impact on the CMB, but they strongly 
suppress the level of mass fluctuations on small scales be- 
cause of the high neutrino momentum before they become 
nonrelativistic. In principle the sensitivity of upcoming sur- 
veys is such tha t it will be possible to test neutri no masses 
below 0.1-leV (Hu, Eisenstein & Tegmark 1998), close to 



those suggested by recen t Super-Kamiokande neutrino re- 
sults (Fukuda et al. 1998). A possible concern is that the ef- 
fect of massive neutrinos becomes important on small scales, 
where the assumption of galaxies tracing dark matter may 
not hold. One of the purposes of this paper is to investigate 
how serious this problem is and, more generally, what is the 



* http://www.mso.anu.edu.au/2dFGRS 
T http://www.astro.prmceton.edu/BBOOK 



relation between the observed redshift space galaxy power 
spectrum and the underlying linear dark matter spectrum. 

The relation between the galaxies and the dark mat- 
ter clustering h as recently been analyzed in the context of 



the halo model (Seljak 2000; Peacock & Smith 200C; Scocci 



marro ct al. 2000). In this model all the mass in the universe 



is divided up into halos of different mass. These halos cluster 
according to the linear theory, up to an overall amplitude 
which depends on the halo mass (halo biasing). To these 
correlations important on large scales one adds correlations 
on small scales, which arise from within the same halos. 
For the latter one needs to specify the radial halo proflle, 
which can also be a function of halo mass. This approach 
has been succesful in reproducing the nonlinear dark mat- 
t er power spectrum and its transition to the linear re gime 
( ^cljak 200d| ; ^coccimarro et al. 2000| ; |Ma fc Fry 2"000| ). 



Galaxies differ from the dark matter in their galaxy mul- 
tiplicity function, which parametrizes the number of galaxies 
inside the halo as a function of halo mass. For a magnitude 
limited sample this function is zero for low mass halos which 
cannot host bright L* galaxies, which already implies that 
galaxies cannot trace dark matter exactly. Above the thresh- 
old the number of galaxies increases with the halo mass, but 
need not grow linearly, as suggested by the gas cooling ar- 
guments where gas in more massive and thus hotter halos 
takes longer to cool and form stars. In addition to the mul- 
tiplicity function there is another effect that changes the 
galaxy clustering properties: one galaxy is expected to form 
at the halo center, which enhances the correlations on small 
scales. These features naturally explain many of the obser- 
vational properties of galaxy clustering in real space, such 
as the power law growth on small scales and the delayed on- 
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set of nonlinear clustering in the translinear regime (Seljak 



200C 



S[;occimarro et al. 2000| ; [Peacock fc Smith 2000^ " 



Previous work in the halo model has focused on the 
clustering in real space, while most of existing and upcoming 
surveys operate in redshift space. Redshift space distortions 
enhance the correlations on large linear scales and suppress 
them on small scales. This changes the clustering pattern 
in a nontrivial way and it is therefore important to include 
these effects when studying the relation between the galaxy 
and the dark matter power spectrum. As an added bonus, 
redshift distortions also allow one to determine f3 — O!^ /b, 
where Q.^ is the mass density and b the bias parameter of the 
galaxies. There are several existing methods to determine 
this parameter in the literature (see ( Strauss & Willick 1991: ) 
for a review). In this paper we address whether different 
methods to determine /3 converge and what is the survey 
size needed for this. 



2 THE HALO MODEL 



where 5c = 1.69 and a{M) is the rms fluctuation in the 
matter density smoothed with a top-hat filter o n a scale 
R = {■AM/A-Kpf''^. We use ( ^heth fc Tormen 199£| ) 

v-1 2p 



b{v) 



and 



1 + 



Scil + v'p) 



(5) 



uf{u)=A{l + u'-^)v'^'^e-''^^ 



(6) 

where p — 0.3 and u' = 0.707i^. The normalization constant 
A is fixed by the requirement that all of the mass lie in a 
given halo 



f{iy)du^l. 



(7) 



On small scales pairs lying within a single halo become 
dominant 



^'^^) ^ W) 



fiu)d. ^ym'- 



(8) 



The halo model uses the Press-Schechter ( Press & Schechtei 



The total power spectrum is the sum of the two contribu- 
tions, 



197'^ ) picture for dark matter, which assumes all the matter 



is in a form of isolated halos with a well defined mass M 
and halo profile p{r, M), which can be modelled as 



Pdm(fc)=Pdm(fc)+Pdm(fc). 



(9) 



p{r) 



Ps 



(r/rs)-«(l-|-r/r^)3+° 



(1) 



where N-body s i mulations give 1 < a < 1.5 (Navarro, Frenk 
fc White 1996|; [Moore et al. 1998|). For the power spec- 



trum analysis adopted here it is convenient to introduce the 
Fourier transform of the halo profile, normalized to unity on 
large scales. 



y{k) 



1 

M 



, , sin(fcr) ,Q 



(2) 



For galaxies the above model needs to be modified in 
several aspects. Instead of the dark matter particles we are 
now counting galaxies inside halos. Small halos cannot host 
very bright galaxies, so there is a lower mass cutoff in the 
halo distribution at a given luminosity cutoff. In addition, 
number of galaxies inside halo need not grow linearly as 
a function of halo mass. Both of these features can be ac- 
counted for by introducing two galaxy multiplicity functions, 
{N){M) and {A''(7V - 1))^''^(M), which count mean number 
of galaxies inside the halo both linearly and pair weighted, 
respectively. Second modification is that radial profile of 
galaxy distribution, t/g(fe), need not be the same as that 



The m ass is determined bv the total mass within the virial of th e dark matter y{k). Both ob servations ([Carlberg et al 



defined to be the radius where the mean density 200C| ) and numerical simulations ( piaferio et al. 1999D show 



radius 

within it is Svir = 200 times the mean density of the universe. 
The concentration parameter c = r^/rs in general depends 
on the halo mass. In this paper we will use a — 1.5 and 
c{M) = 6{M/M^,)'°-'^^, where M* is the nonlinear mass 
scale defined belo w. This choice fits well the results of N- 
body simulations ( Moore et al. 1998 ) and has been shown 
to give good agreem ent with rea l space power spectra from 
N-body simulations ( Seljak 200C ) , but we note that other fits 
with a = 1 and a different choice of c{M) can give equally 
good agreement with these. 

In the halo model the power spectrum consists of two 
terms. The first is that due to a system of correlated halos, 
with inter-halo correlations assumed to be a biased sam- 
pling o f Pti„(k). Since the real space convolution is simplv a 



that at least for some types of galaxies the two functions 
must differ. Here we will for the most part adopt the ap- 
proach where the two are equal, but we also discuss the 
modifications when this assumption is dropped. Finally, one 
expects there will be one galaxy which forms at the cen- 
ter of the halo. The correlations between this galaxy and 
the rest of the galaxies inside the halo will be sensitive only 
to a single convolution in the radial profile. For large halos 
with {N{N — 1))^ ^ 1 the presence of the central galaxy 
does not change significantly the number of pairs or their 
statistics. For small halos where {N{N — 1))^ ^ 1 its ex- 
istence changes the correlations significantly and is in fact 
necessary to explain the steep pow er law in the galaxy cor- 



relation fun ction to small scales (Seljak 2000; [Peacock & 



Fourier [space multiplication this contribution is 



Smith 2000) 



ega 



P''\k) = Pii„(fe) 



f{u)dub{u)y{k;M) 



Putting the above together we have. 



(3) 



P,,{k) = PUk) 



where b{u) is the (linear) bias of a halo of mass Miy) and 
fiy) is the multiplicity function. The peak height v is related 
to the mass of the halo through 



(4) 



I j fi^)d^^-^b(i.)y,(k,M) 



a{M) 



where n is the mean density of galaxies in the sample, 

-W'f^''^'''' = -p- 

The Poisson term is given by 



(10) 



(11) 
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PL{k) 



(27r)3n2 






1)> 



y{k, Mr, (12) 



where we approximate the effect of the central galaxy by 
using p — 1 for {N{N — 1)) < 1 and p = 2 otherwise. 

In redshift space there are two effects which modify the 
above expressions. The first is a boost of power on large 
s cales due to streaming of matter into overdense regions 
( Kaiser 1987 ) . The second is a reduction of pow er on small 
scales due t o virial motions within an object (Peacock & 
1994). In the halo model these two effects can 



Dodds 



be 



separated into the halo-halo and one halo (Poisson) contri- 
butions. In this model the virial motion suppression becomes 
a function of scale, since larger halos (with larger velocity 
dispersions) dominate at larger scal es than smaller halos . 
In contrast to the previous models (Hatton & Cole 199!:) 



t his model su ccesfully reproduces N-body simulation results 
( [White 2000| ). 

In linear theory a density perturbation Sk generates a 
velocity perturbation S — —ikv with v parallel to k. Using 
the plane-parallel approximation (fcr S> 1) the redshift space 
galaxy density perturbation 5^^ is given by 



^r = ^g + ^vA* I 



(13) 



where /j. — f ■ k, Sg is the real space galaxy density per- 
turbation and 5v is the velocity divergence. This can be re- 
lated to the density perturbation 5dm via 5^ = /5dm, where 
f{fl) = dlogS/dloga ~ Q,^ and a is the scale-factor. 

On small scales virial motions within collapsed objects 
act as a gaussian convolution in redshift space, which sup- 
presses power. We will model this as a gaussian filter with 
mass dependent 1-d velocity dispersion a, acting on the 
mode component along the line of sight. Assuming that the 
halos are isothermal we may use the mass within the virial 
radius to obtain the ID velocity dispersion of a halo of mass 
M, 



a = [GM/2r, 



il/2 



THor^iT, 



(14) 



where Ho is the Hubble constant and r^ir is the virial ra- 
dius, which can be related to the halo virial mass, Mvir ~ 
47r(5virp7'5ir/3.The density contrast in redshift space is 



5„ = S^e 



'(ka-iJ.)^/2 



(15) 



Sirlce the power spectrum is in general a functinn of 
both k and /x we must decide which quantity we are inter- 
ested in before proceeding. The most common is to average 
over n to obtain the isotropized power spectrum. On large 
scales this implies averaging over ^ the square of equation 
hA. This can be further improved by including the small scale 
dispersion. A choice that seems to work well in comparison 
to simulations is to add a term obtained by radially averag- 
ing equation ha (White 2000). On small scales we only need 
to average over the square of equation [l5l since there is no 
linear effect. 

Combining the above the isotropized redshift space 
power spectrum in the halo model becomes 



^o(fc) = (Fg + |i^vfg + ^F^') Plin(fc) 

{iV(iV-l)) 



-f 



(27r)3f7.2 



— f{v)du- 

p ' 



KP 



(16) 
7^J,(fca)|^/g(fc,Af)|^ 



where 



Fv = f f{u)dub{iy)TZi{ka)y{k;M) 

^g = f f f{u)dJ-^^^b{u)n,{ka)ygik,M) (17) 



and 



Tlr, 



ka[p/2] 



1/2n 



y'vF erf(a) 



(18) 



for p = 1, 2. 

Figure pi shows various bias functions as a function of 
k, defined as the square root of the ratio of the galaxy to 
the linear dark matter power spectrum. Throughout we use 
A.CDM model with Q,m = 0.3, Q,x = 0.7, normalized to 
erg = 0.9 today. Linear dark matter power spectrum is used 
here, since this is the quantity we wish to reproduce. Also 
shown is the square root ratio of the nonlinear redshift space 
dark matter spectrum to the linear power spectrum. This 
quantity can be compared directly to the N-body simula- 
tions. We verified it to be in a remarkable agreement with 
those (see also White 2000). On large scales it should cor- 
respond to 6 = 1.18. It is interesting to note that the red- 
shift space dark matter power spectrum agrees to within 
20% with the real space linear power spectrum over the 
entire range of scale. This is somewhat coincidental, since 
for the nonlinear redshift space power spectrum the corre- 
lated halo-halo term is suppressed on small scales because 
of virial motions (TZ\{ka) < 1) and because of halo profile 
{y(k;M) < 1). The difference is picked up by the Poisson 
term, which does not depend on the linear power spectrum, 
except through the mass dependence of the concentration 
parameter. 

For galaxies we see that the linear bias, defined as the 
ratio of the redshift space galaxy power spectrum to the 
linear dark matter spectrum, typically exceeds the redshift 
space dark matter bias on large scales. This is of course not 
surprising and reflects the fact that galaxies are a biased 
tracer of dark matter. For galaxies found predominantly in 
groups and clusters the bias is larger than for those which are 
also found in the field. The choice of the galaxy multiplicity 
function used here is mo tivated by the sem i- analytic mod- 
els of ga laxy formation (Benson et al. 2000|; [Kauffmann et 



al. 1999). Least biased are regular galaxies selected only on 



the basis of their luminosity (dashed line). More biased are 
red galaxies with Mb — My > 0.8 in semi-analytic models 
of Kauffmann et al. (1999; dash-dotted), while most biased 
are those with Mb — My > 2 (dotted). It is important to 
emphasize that these are just plausible choices of the galaxy 
multiplicity function as the data at present do not allow one 
to determine these directly. Theoretical models of galaxy 
multiplicity function can vary at least at the level of 30%, 
leading to variations of up to 150 km/s in 1-d velocity dis- 
persion. For this reason we only emphasize the features of 
the model which are generic and expected for any choice of 
the parameters, even though the relative importance of dif- 
ferent effects may vary from model to model. For example, 
in all cases the galaxy multiplicity function is expected to in- 
crease less rapidly with the halo mass than the mass itself, 
but may have a different low mass cutoff and/or different 
shape. 
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The bias is scale independent on large scales for the reg- 
ular galaxies. For these bias becomes scale dependent above 
k > O.l/iMpc"^, on scale somewhat smaller than for the 
dark matter. This is because these galaxies are preferentially 
found in smaller halos relative to the dark matter, and are 
hence dominated by systems with smaller a, thus the virial 
motion suppression is smaller. Just as in the case of the dark 
matter the galaxy linear bias first declines with scale. This 
is because correlations are suppressed both by the finite ex- 
tent of the halos and by the virial motions within them. At 
even smaller scales, above k > l/iMpc~^, bias begins to rise 
again. This is caused by the nonlinear Poisson term for the 
galaxies, not present in the linear power spectrum, which 
enhances the correlations on small scales. This is of course 
the term that gives rise to the nonlinear clustering pattern 
and the familiar power law slope. ,]ust as in the real space in 
the redshift space the enhancement is also more important 
for the galaxies than for the dark matter. One reason for this 
is that galaxies are preferentially found in smaller systems 
relative to the dark matter, so the suppression because of 
finite halo extent and virial motions is less important than 
for the dark matter. Another is that on smaller scales the 
correlations are dominated by small halos with (A'^) < 1. In 
these there is a central galaxy, which does not contribute to 
the suppression because of halo profile and virial motions, 
leading to an enhancement of the galaxy power spectrum 
over that of the dark matter. Finally, a more concentrated 
distribution of galaxies relativ e to the dark matter , as sug- 
gested by some observations (Carlberg et al. 200C), would 
also lead to an enhancement of the galaxy power spectrum 
relative to that of the dark matter. 

The choices motiv ated by red or elliptical galaxies in 
semi-analytic models (Kauffmann et al. 1999) elliminate 

Salaxies in less massive halos. These are shown in figure 
as the dash-dotted and dotted lines, for the case of red 
{Mb - Mv > 0.8) and very red [Mb - My > 2) galaxies, 
respectively. The bias in this case shows scale dependence 
even on very large scales, k < O.Ol/iMpc"'^. This is because 
the Poisson term becomes important again relative to the 
linear power spectrum on very large scales, where the latter 
approaches slope n ~ 1, while the former remains at n = 0. 
Galaxies that are only found in rarer, more massive, systems 
have the Poisson term that is larger, so this effect is rela- 
tively more important for this type of galaxies as opposed 
to the normal galaxies. As pointed out in Seljak (2000) this 
effect is not present for dark matter because of mass and mo- 
mentum conservation, which requires the power spectrum of 
a local process to decline as P{k) oc fc* on large scales. Since 
the Poisson term is constant on large scales one can attempt 
to model it as such and remove it. Note that its amplitude 
is larger than the usual Poisson term arising from the dis- 
cretness of galaxies, which depends on the total number of 
galaxies. Here the amplitude of Poisson term is determined 
by the effective number of halos, which is smaller than the 
number of galaxies. If one removes this constant term then 
the remaining power spectrum should trace dark matter, so 
bias should become scale independent again. 

For the same reason that these galaxies are in more mas- 
sive halos the virial motion and finite halo size suppression 
become important on larger scales and the bias is already 
scale dependent for k ~ 0.1/iMpc~^. On smaller scales the 
bias dependence is also sensitive to whether the galaxies are 




0.001 



Figure 1. Bias b{k) defined as the square root of the ratio of 
redshift space nonUnear power spectrum to tlie linear real space 
dark matter spectrum. From top to bottom are shown very red 
galaxies {Mb — My > 2, dotted), red galaxies {Mb — My > 
0.8, long dash-sotted), normal galaxies (short dashed) and dark 
matter (solid). The bottom two lines are ratios between massive 
and massless neutrino transfer functions, for rrii, = O.leV (upper, 
short dash-dotted) and m^ = leV (lower, long-dashed). 



central inside the halo or not and whether the radial distri- 
bution follows that of dark matter, both of which can have 
important impact on the bias. For example, if red galaxies 
avoid centers of halos (because the gas preferentially cools 
to the halo center, where it can form new stars which are 
blue), then the bias is a less rapidly rising function of k on 
small scales than that shown in figure nl) . 



3 /3 FROM REDSHIFT DISTORTIONS 

Another application of the redshift distortion analysis is to 
extract the parameter (3 — f/b. There are many different 
ways to do this (see Strauss & Willick 1995 for a review). 
The two examples used here are by combining the galaxy 
power spectrum Pgg with the velocity power spectrum Pvv or 
their cross-spectrum Pvg and by using Legendre expansion of 
redshift space power spectrum. If bias is constant on large 
scales /3 will also be a constant, but in the more general 
case considered here it will be scale dependent. Moreover, 
different methods may not even agree on a given scale, so 
the meaning of l3 itself becomes questionable. 

Given that in general the power spectrum is a function 
of both fi and k it is not possible to extract the three com- 
ponents Pgg, Pvg and Pvv from it uniquely, since it depends 
on the adopted procedure. In fact, a clean separation into 
the three components is only possible on large scales where 
only the linear compression redshift distortion operates. To 
obtain an idea what happens when small scale effects be- 
come important we perform angular averages assuming lin- 
ear theory model and reconstruct the 3 components from 
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these. As before we assume plane-parallel approximation. 
For isotropic averaging in equation [17| this gives 

Po(fc) = i y" d^^ {S,+ ^j.^5,Y = Pgg + |Pvg + iPvv.(19) 

Similarly we can also perform averaging with fj,^ and fj."^ 
weights, 

1 /•+! 9 1 9 1 



1 /"^^ 1 






From Po, P2 and P4 we can reconstruct uniquely Pgg, Pvg 
and Pvv using the above expressions. At the same time the 
halo model also gives predictions for the averaged spectra 
P2 and Pi just like for Po (equation [17^). These are given by 

P^ik) = [^Fl + ^F.F, + ip,^) PUk) 
M {N{N-1)) 



+ 



(27r)3n2 



— f{v)dv^ 



KP 



TZp+2{ka)\y^{k,M)f 



P^k) = (iPg + yPvFg + iPv') PUk) 
Af ,, , , {iV(Af-l)) 



(27r)3n2 



— f{y)du- 



M2 



where 



TZp+q{ka) — / /i''e 
'o 



<j^-p(fc<7M)V2 



dfj., 



1lp+4{ka)\y^{k,M)\'\ 
(21) 



(22) 



which h ave simple analytic e xpr e ssiong similar to equation 



to Pgg, Pvg and Pvv Their reconstructed scale dependence 
of /3 show a similar behaviour as our model, where /3 from 
Pvg/Pgg declines with k, while that from (Pw/Pgg)^ in- 
creases with k (see their figure 4). We caution that this com- 
parison is just illustative, since the two methods of analysis 
differ in details (such as the use of plane-parallel approxi- 
mation) and so cannot be directly compared. 

We can also extract the cross-correlation coefficient 
r = Pvg/[PvvPgg]^ from this analysis. For normal galaxies 
it is close to unity for k < 0.1/iMpc^^ and rapidly declines 
above that. This means that the galaxy density and veloc- 
ity divergence as reconstructed from this method become 
poorly correlated on small scales. Given this it is meaning- 
less to combine the different estimates of /3 to enhance the 
statistical significance, since they do not measure the same 
parameter and the two fields are only poorly correlated on 
small scales. 

For biased (red or elliptical) galaxies found predomi- 
nantly in groups or clusters the correlation is smaller than 
unity even on large scales and the agreement between differ- 
ent reconstructed /3 is less good there. This is mostly caused 
by the Poisson term, which is not negligible for galaxies on 
large scales. Divergence between /? from different methods 
appears already for k ~ 0.04/iMpc~'^, caused by the more 
significant influence of massive halos with larger velocity 
dipersions. One should therefore be specially careful when 
drawing conclusions on cosmological parameters from this 
sample of galaxies, as for example from the planned bright 
red galaxy sample (BRG) in SDSS. 

Another often used way to extract /3 is to determin e 
the ratio of quadrupole to monopole terms (Kaiser 1987), 
which in terms of the above quantities is given by Pgm = 
2.5(3P2/Po — 1). In linear theory one predicts it to be ( |Cole 



Fischer fc Weinberg 1994) 



7^p+2(a = ka[p/2Y'') = ^ 



TZp+4{a = ka[p/2Y^') = 



30F 




for p = 1, 2. 

The procedure to extract /3 is the following: first we 
compute Po, P2 and P4 from above expressions. Next we as- 
sume they are determined by the linear combinations of Pgg, 
Pvg and Pvv as valid in linear theory, which allows one to de- 
termine them uniquely. We then take ratios [Pw/PggY and 
Pvg/Pgg to determine /3. These results are shown in figure 
for normal (top) and red (Ms — My > 0.8; bottom) galaxy 
sample used in figure fD. We see that for normal galaxies 
the two reconstructed (3 functions are approximately con- 
stant and equal for k < O.lftMpc"^. For k > O.lftMpc"^ the 
two /3 functions diverge away from the large scale value and 
away from each other. This indicates that on scales below 
50/i~^Mpc one cannot extract the true value of /3 and that 
different methods of determining it can give rather different 
answers. This is because on these scales virial motions within 
halos become important and pure linear theory ansatz is no 
longer valid. The scale dependence of the two (3 reconstruc- 
tions is in qualitatively good ag reement with the behaviour 



_ 1 + 2/3/3 + /?V5 
■""" 4/3/3 + 4/32/7 ■ ^^^> 

By solving the quadratic equation above one can determine 
/3 as a function of scale. This quantity is also shown in fig- 
ure y and has a similar behaviour to /3 from Pvg/Pgg ratio, 
although the suppression on small scales is delayed relative 
to it. 



4 CONCLUSIONS 

We have analyzed the redshift distortion effects using the 
halo model, which has proven to be remarkably succesful 
in explaining nonlinear re al space power spectrum of both 
galaxies and dark matter ([Via fc Fry 200C ; Peacock fc Smith 



200c ; Scoccimarro et al. 2000; seljak 2000), as well as 



redshift space power spectrum of the dark matter ( White 
200c). We have shown that redshift space bias for regular 
galaxies is likely to be constant only for k < O.l/iMpc"'^. 
Between 0.1/iMpc~^ < k < l/iMpc~^ the bias declines by 
10-20% and rises again above that. This may be important 
for attempts to extract th e value of cosmological parameters 



from such measurements (Tegmark, Zaldarriaga & Hamilton 



seen in real data . For example, (Hamilton, Tegmark & Pad- 
manablian 2000|) have decomposed the data in a similar way 



200c ) , which typically assume nonlinear effects are negligible 
up to fc ~ 0.2 - 0.3/iMpc"\ 

One parameter that is particularly sensitive to this ef- 
fect is neutrino mass, which suppresses the power spectrum 
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correlation coefficient. Top figure shows regular galaxies found 
also in the field, bottom shows biased (red, elliptical etc.) galax- 
ies found predominantly in groups and clusters. 



roughly in the same scale range where bias also becomes 
scale dependent. Since it has the same effect the two can 
be degenerate over the quasi-linear regime and this would 
complicate the attempts of accurately determining the neu- 
trino mass from such galaxy clustering measurements. Fig- 
ure hi shows that the effect of neutrino mass begins to affect 
the power spectrum already on scales larger than where the 
bias is scale dependent. This is good news for the efforts 
to extract the neutrino mass, but the effect can only be 
measured if very large scale correlations can be accurately 
measured. These are limited by finite volume sampling vari- 
ance, which can only be reduced by having a larger vol- 
ume, so that the amplitude of correlations can be reliably 
determined for k <^ O.l/iMpc"'^. If only the information 
around k ~ 0.1/iMpc~^ is used then one cannot separate 
the scale dependent bias effects from those of massive neu- 
trinos. Current surveys such as PSCz still have large errors 
for k ^ 0.1/iMpc~^ and the new generation of surveys, such 
as SDSS and 2dF, combined with a more detailed modelling, 
will be needed to determine the neutrino mass from these 
measurements. 

Galaxies found in rare systems such as groups and clus- 
ters suffer from another effect. For such galaxies the Poisson 
term can be so strong that it can exceed the linear corre- 
lation term not only on small scales, but also on very large 
scales. This effect is not present for dark matter, which is 
protected fro m it by causalit y and conservation of mass and 
momentum (Zeldovich 1970). Because of this the bias can 
rise again for k < 0.01/iMpc^^. Since on such large scales 
this is a pure Poisson term with the slope exactly n — 
one can attempt to model it as a sum of two contributions 
and its amplitude can be estimated from the small scales. 



Equivalently, one can perform the correlation function anal- 
ysis, where this term is not present on large scales. 

Determination of parameter /3 shows similarly a con- 
veregence to a single value for k < 0.1/iMpc~^, at least 
for normal galaxies where the Poisson term does not be- 
come important on very large scales. Around and above 
k ~ 0.1/iMpc~^ /9 rapidly becomes scale dependent. The 
actual behaviour depends on the specific analysis and /3 can 
either grow or decline with scale. This is caused by the effect 
of virial motions, which counter the linear compression effect 
and are more important in more massive halos, which dom- 
inate the nonlinear clustering on large scales. The model re- 
produces well the scale dependence of /? seen in the analysis 
of the real data (Hamilton et al. 2000) . The cross-correlation 
coefficient between the galaxy and the velocity field diver- 
gence as obtained from the redshift space distortions is close 
to unity below k ~ 0.1/iMpc~^ and rapidly declines above 
that. Galaxies and velocities are singificantly less well cor- 
related than the dark matter and the galaxies in real space 
analysis, where the cross-c orrelation co efficient is unity at 
least up to fc ~ l/iMpc~^ (Seljak 2000). Our results are in 



broad agreement with other recent analysis of non linear bias 
and its eff e ct on redshift space distortions and /? (Hatton & 



Cole 199c 



Berlind, Narayanan & Weinberg 2000). Together 



these results argue for a need of a more refined analysis 
of redshift distortions and redshift space power spectrum if 
the statistical power of existing and upcoming redshift space 
surveys is to be fully exploited. This is needed both to reduce 
the systematic effects and to extend the analysis to smaller 
scales. Since the halo approach used here reproduces well 
the results of N-body simulations and semi-analytic models 
it can serve as a useful framework within which one can ex- 
tract the true underlying cosmological model from the data. 
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